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Abstract 

We derive semiclassical quantization for a spin, study it for not too small a spin quantum number (S > 
5), and compute the 25+1 eigenvalues of a Hamiltonian exhibiting resonant tunneling as the magnetic 
field parallel to the anisotropy axis is increased. Special attention is paid to the resonance condition. 
As a corollary we prove that semiclassical quantization and quantum-mechanical perturbation theory 
agree there where they should. 
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1 Introduction 



Resonance has revived the experimental interest in spin quantum tunneling. The idea is both fascinating 
and simple Each molecule in, e.g., a Mni2 acetate crystal carries a spin of fixed angular momentum 

f Oh and experiences a constant magnetic field H. The corresponding Hamiltonian is given by 

H = -jS 2 - 5M bH • S . (1) 

The first term on the right is a magnetic anisotropy. Now let H = (a, 0, d) and absorb gfiB into a and 8 
so that Eq. (|l|) reduces to 

U = -[jS 2 + SS Z ] - aS x = -F(S Z ) - aS x . (2) 

Here 7 > and we may assume a, 5 > 0. By putting 7 := —7 we encounter a situation whose tunneling 
physics and mathematics hardly changes and therefore need not be treated separately. In the present 
paper, S x and S z have the dimension of angular momentum, to be measured in units h. For the moment 
we suppose that the spin quantum number S is an integer. 

The term aS x tries to generate a rotation about the cc-axis and thus aims at inducing a tunneling 
transition. It certainly does so for 6 = 0. For arbitrary nonzero 8 the degeneracy in the spectrum of S 2 
is lifted and no tunneling can occur unless a special choice of S is made that restores the degeneracy. For 
a = a degeneracy exists, if F(mh) = F(nh) for some integers m and n, i.e., putting 7?L = T , 

fiW = -m 2 T - mS = ~n 2 T -nS= . (3) 
In passing we note that both V and 6 have the dimension of frequency. Equation (^) tells us 

(m 2 - n 2 )T + (m - n)5 = (4) 
so that with — S < m n < S we are left with 

(m + n)r + S = =>■ m + n = -6/T = -k (5) 

where k > 1 is chosen to be a positive integer. Hence m = —n — k. The remarkable aspect of this condition 
is that it can be realized for all n in the ballpark —S<n< S—k. Furthermore, degeneracy can be realized 
by a single choice, viz., i5 = kT. Here one has exploited the fact that the anisotropy is quadratic in S z - 
The variable parameter 8 is at the disposal of the experimentalist and the corresponding magnetic field 
is usually increased from minus to plus a few Tesla. One then finds several resonances in between. 

It is fair to say that the above setup is indeed elegant. It is a corollary of semiclassical quantization 
as treated in the present paper that the resonance condition (|^) also holds for a ^ 0. The (tiny) level 
splitting is a consequence of quantum mechanics. 

Meanwhile resonance has been studied intensively. Already those who invented it initiated a first 
attempt [0 at its theoretical understanding. In an admirable tour de force, Chudnovsky and Garanin 
H were able to fully analyze the influence of the heat bath provided by the surroundings of the spins 
on the tunneling process if <5 = 0, and there is little doubt that their arguments can be generalized to 
the 5^0 case. Here we assume a much lower temperature so that the heat bath can be neglected. The 
analysis of semiclassical quantization, including resonance, is the main theme of the present paper. It 
yields the energy eigenvalues of the Hamiltonian (^) to a good approximation and thus completes our 
earlier work |j, ^ on the WKB formalism for spins; we refer the reader to the excellent review pi 
for additional information regarding the experimental and theoretical context. In §2 we reconsider the 
WKB formalism and indicate the use of the WKB wave functions in computing the level splitting at 
resonance. We then treat semiclassical spin quantization (§3) and determine (§4) the energy eigenvalues 
for arbitrary external field, i.e., 6. In §5 we extend the results to half-integer spin quantum numbers and 
Hamiltonians in which the term aS x is replaced by aS l x with Z > 1. The conclusion is that semiclassical 
quantization holds here too but that for even I the condition (|^) may correspond to true degeneracy; that 
is, no resonance and no tunneling. We finish the paper with a discussion (§6). 
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2 WKB and Resonance 



The eigenvalues and eigenvectors of the Hamiltonian (|^) are solutions of the eigenvalue equation Tltp = 
Eip. The S z axis being the main anisotropy axis, it is natural to write this equation in a representa- 
tion with S z diagonal. Then S z is simply s — nh, a multiplication operator on the spectrum of S z . 
Furthermore, let T±^ induce a translation by ±h so that (T±nip)(s) — ip( s i h). Then S x reads 

S x = \ (S+ + S-) = ^ [a + T h + a (^(s-fi)) T_ t ] (6) 

where a(s) — [a {a + K) — s 2 ] 1 / 2 . As a consequence, the Schrodinger equation assumes the form of a 
second-order difference equation, 

h n , n -xi>n-l + h n ,n+11pn+l + (h nn ~ E)ljj n — . (7) 

The matrix elements and the vector components are taken in the basis of the eigenvectors of S z , and (Q) 
is valid for — S+l<n<S — 1. Disregarding the remaining two equations, relation (|7|) has two linearly 
independent solutions for any value of E; they are determined by fixing, e.g., ipo and ip^. The 2S + 1 
eigenvalues of Ti. are singled out by requiring that tp satisfy the boundary conditions 

h±s,±(S-i)^±(s-i) + (h±s,±s ~ E)ip± s = . (8) 

Equivalently, we can extend (Q) to n = ±S by defining arbitrary real /i±s,±(S+i) = ^±(s+i).±s an( i 
imposing the boundary conditions ip±(s+i) = 0. 

If a = 0, the eigenvectors are those of S z and the eigenvalues equal with n ranging from 

-S to S; cf. Eq. (|). It is interesting to note that, as the energy increases while avoiding accidental 
degeneracies (||), the eigenvectors are alternately localized on the 'right' (n > 0) or on the 'left' (n < 0). 
Far from degeneracy an approximate localization on alternating sides remains true for a nonvanishing a, 
if a <C TS + 5, which is supposed throughout the paper. 

By varying d, neighbouring eigenvalues can get very close to each other, although true degeneracy 
cannot occur |l2]. As we discuss below, the mechanism of the avoided level crossing is spin tunneling. 
When two eigenvalues become as close as possible, we speak of a quantum-mechanical resonance. This has 
to be distinguished from what we call a semiclassical resonance, which we understand to be the coincidence 
of semiclassical eigenvalues. We shall prove in §4 that the condition for a semiclassical resonance remains 
(||). Apart from the resonance at 5 = 0, we cannot expect that the two definitions predict the same 
'resonant' values for S. Because, however, semiclassical eigenvalues nicely approximate the true ones, 
only a small 'fine tuning' of 6 may be necessary to pass from semiclassical to true, quantum-mechanical, 
resonance. In the remaining part of this section we explain how one can estimate the level splitting at 
resonance. 

Let us choose 5 close to a resonant value and let E\ < E% be two (unknown) neighbouring, nearly 
degenerate, eigenvalues. Let us also imagine that we are given two linearly independent vectors r\(E) n 
and •d{E) n which depend continuously on E and solve (0). To avoid all confusion, we emphasize that 
no linear combination of them satisfies the boundary conditions (||) and, hence, the eigenvalue equation, 
if E is not an eigenvalue. On the other hand, the eigenvector belonging to E\ is a linear combination 
of rj(Ei) and $(Ei), and its analog holds for E%. Let e[ sc ^ and E^ be the semiclassical eigenvalues 
corresponding to E\ and E2, respectively. Because all four energies are now close to each other, we can 
replace r](Ei) and 77(^2) by n^E^), and 'd(Ei) and ^(£'2) by ^(E^), or vice versa. In this way, finding 
Ei and E2 and the corresponding eigenvectors reduces to good approximation to diagonalizing H in a 
two-dimensional subspace. In so doing, we can use the WKB method to obtain n{E) and "0(E). 

Though Wentzel, Kramers, and Brillouin (WKB) never considered spins, their idea is applicable here 
too, provided one generalizes the formalism appropriately (8[ |9[ [l0| so as to take care of the discrete 
nature of a spin and its different commutation relations as compared to a particle. A semiclassical 
analysis formally means that we take the limit fi — > and at the same time S — > 00 in such a way that 
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Figure 1: Classical orbits of the Hamiltonian H = — 7>Sf — 5S Z — aS x are the intersection(s) of the energy 
surface E = —~fS z — 6S Z — aS x and the sphere 5j + Sy + S z = a 2 , here plotted for S — 20, ft = a = 7 = 
1, 5= 10, and E = -50. 



HS = a remains constant. In this limit we are left with a continuum description with s ranging in the 
interval [— cr, cr]. In the spirit of WKB we now make the ansatz ip — exp(iS/H) with 

s = s + y2(-rs n (9) 

n=l 

for the wave function we are looking for, expand everything in powers of ft, and usually stop after the 
first-order term, the zeroth-order one being dominant. So de facto we use a continuum description, even 
though k is still finite. 

The dominant contribution to tunneling comes from the classically forbidden region, say, between the 
inner turning points b\ < 62- In this interval the two linearly independent WKB wave functions read 
110 

1 f s ( —E — Fis'W 

4>e,i{s) = C;exp-~ / ds'arccosh ( — — I (10) 

ftj bl \ aa(s') J 

and 

1 f b2 f-E- F(s')\ 

4>e r{s) = C r exp — — / ds'arccosh — - — . (11) 

nJs V aa{s') J 

Here C/ ir are normalization constants and the arccosh expression stems from 1S0 in (^). As can be seen 
from Figs. 1 and 2, there are four turning points on the s :— S z axis, two inner ones, b\ and &2j and 
two outer ones, a\ and 0,2- The inner turning points b\ and 62 are boundaries of the classically allowed 
motion to the left and right of them, respectively, and such that here the argument of the hyperbolic 
cosine equals 1. The functions depend continuously on the (classical) energy E. By using appropriate 
connection formulae, both can be extended to the whole interval [— cr, cr], and the functions obtained in 
this way are almost perfectly localized on the left and on the right, respectively. Indeed, between b\ and 62 
Eqs. (10) and ( JIl| ) define, respectively, a rapidly decaying and a rapidly increasing function. Localization 



of the extended solutions then follows from 4>E,i{bi) ^> 4>E,l(b2) and 4>E,r(p\) <C 4 > E,rQ)'z)- F° r later use 
we note that the two functions are not orthogonal to each other and, because of the localization, their 
tiny overlap comes essentially from the classically forbidden region. 

If we restrict 4>e,i(s) and 4>E, r {s) to the discrete values s = nh with integer n between —S and S, 
we obtain the two vectors rj{E) and §(E). Our earlier discussion shows that, if E happens to be an 
eigenvalue far enough from other eigenenergies, then 4>e,i or 4>E,r alone is a good approximation of the 
corresponding eigenvector. Returning to the problem of resonance, we can find E\ and E2 by diagonalizing 
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Figure 2: The effective, quartic double-well potential U(s) given by Eq. ( |2l| ) has been plotted for 
% = a = 7 = 1, 8 = 10, and E = -56.7. For S = 20, E = E_ u = E A ; that is, the levels -14 and 4 are at 
resonance. The two horizontal segments are located at the 'energy' level £ = (S 2 — E 2 )/2. They end up 
in the turning points, which limit the classical motion. U" is about 1% bigger in the minimum on the 
right than in the minimum on the left. 
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TL in the subspace spanned by, say, <p E {^) ; and </> E (sc) r . Close to the semiclassical resonance (|) there 
exist integers m < < n such that these functions are localized near mh and nh, respectively. We 
shall use the shorthand (f> m and 4> n for them. Let x denote the matrix of TL so that Xij = { < Pi\7~L\4>j) i 
and let o be the overlap matrix with elements Oy ~ (</>i|</>j), where i,j 6 {ro, n}. We recall that o mn is 
nonvanishing, although very small; we shall estimate it below. There is no harm in supposing <j) m and <p n 
to be normalized. 

We now turn to an explicit calculation of the quantum-mechanical level splitting. To a good approx- 
imation ||, ||, [To) , Ei and E2 agree with the eigenvalues E- and E + of the 2x2 matrix 

q — 1 J Xmm XmnOmn \mn X^n^m/i j (12) 

1 — <^mn V Xmn Xmrn^mn Xnn XrnnOmn J 

To see it, expand the eigenfunctions in terms of 4> m and 4> n and take matrix elements. For the eigenvalues 
we find ^ 

E± = — 5 — AXmm + Xnn-1XmnOmn±D) (13) 

2(1 - 2 mn ) 

where 

1/2 

£>= [(Xmm-Xn«) 2 +4(X )(x (14) 

Let (a;±, denote the eigenvectors of o _1 x corresponding to That is, the approximate eigenvectors 
of TL are x±<p m + y±4> n . For r± — x±/y± an elementary computation yields 

„2 Xn " " S± (15) 



Xmm -E± 

and, by orthogonality, r + r_ = — 1 — (r + + r_)o TOn . Near resonance Eqs. ( [l3"l ) and ( |l5| ) provide improved 
eigenvalues and eigenvectors as compared to the semiclassical ones. 

We now have to tune 8 to quantum-mechanical resonance. Supposing Xmn is of the order of the 
overlap, the second term under the square root in ( jli] ) is of order o 2 mn . The overlap is a smooth function 
of 8 and remains uniformly small in the small region where we vary 8. So the minimal distance between 
E+ and £L is reached when A mn — Xmm — Xnn vanishes (or is also of order o mn ). Thus we conclude that 
the quantum- mechanical level splitting is of order o mn . The order of magnitude of o mn is easily inferred 
from Eqs. jl0| [Tl| ). For this estimate we may suppose E = e[ sc ^ = E^- Then 

Y r^2 / F(s)\ 

Omn = Cexp — — / ds arccosh ( ) (16) 

hj bl V aa ( s > J 

where C is a constant of order 1. The exponential factor in o mn can be interpreted as a transition 
probability. This suggests to write the level splitting AE = E% — E\ in the form 

. _ tt/& 1 f b2 , , f-E-F(s)\ 

AE = — exp-- / ds arccosh (17) 

r hj bl \ aa{s) J 

where 1/tq is an attempt frequency. Recalling the expression of the hyperbolic cosine in terms of a 
natural logarithm and using b\ pz mh and 62 ~ nh, we find the order-of-magnitude estimate 



. „ trh 
AE « — 

to 



aS 



r(m 2 + n 2 ) + 8{m + n) 



(18) 



By identifying tq with the time period of the classical motion - see Eq. (53) below -, Eqs. (|l7|) and ( |l8| ) 
become fully explicit. For 8 = and n = — m = S* they agree, respectively, with formulas (C.ll) and 
(C.12) in Appendix C of our earlier wor k l| . In that case dl^) leads to a remarkably precise result, as 
can be seen in Table 1 of ||. Equation (p.7[) was obtained ]lCfl also by an independent argument. 

Of course 8 can be chosen to let A mn vanish. The reason is that Xmm ~ E[ bC ^ and Xnn ~ E^ (the 
small deviation coming from the fact that <\> m and <f> n do not satisfy the boundary conditions (|8j)), and 
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semidassical eigenvalues do cross each other at resonance. Thus A mn — for a S close to — (to + n)T. 
In particular, because of a reflexion symmetry of the Hamiltonian, the level pairs {m = — n, n} are at 
semiclassical and quantum-mechanical resonance, once S — 0. 

The origin of the level splitting is spin tunneling. At quantum-mechanical resonance \r±\ — 1 and 
thus the approximate eigenfunctions are ifi± = (<f> m ± <p n )/^/2. If we start on the left, we take (f> m = 
{tp + +V-'-)/v / 2. It evolves under the influence of the dynamical evolution generated by exp(itH / h) . After 
a time T given by Eq. (17) through TAE — nth, the system is in <f> n = — ^>_)/\/2, i.e., on the right. 
All this is exactly as in the case with a reflection symmetry S z <-> — S z such as when S = 0. In the 
corresponding classical problem an approximate reflection symmetry survives for remarkably high values 
of S (5 < 10r) with a shifted center of symmetry; cf. Figure 2. Because of this approximate symmetry 
and for low enough energies, the attempt frequencies in the two, now different, orbits centered at to and 
n are hardly different; see Eq. ( p2| ) below. This fact is crucial for the interpretation of tq in the level 
splitting formula (|l7|). 



3 Semiclassical Quantization 



Semiclassical quantization of a single spin can be handled straightforwardly since we always find closed 
orbits - if any - as the intersection of the energy surface Ti = — [yS z + SS Z ] — aS x = E and the sphere 
S 2 = S 2 + Sy + S 2 = er 2 ; cf. Fig. 1. Here the classical equations of motion associated with the Hamiltonian 
H = —F(S Z ) — aS x can be reduced to a second-order differential equation O] for S z , 



S Z = -F(S Z )F'(S Z )-EF'(S Z ) 



! S Z = 



U(s) 



(19) 



S=S; 



where 



In the present case, F(s) = 7s 2 



U(s 
Ss so that 



If 2 (s)+EF( S ) 



1 2 2 

-as 
2 



(20) 



U(s) = i 7 2 s 4 + jSs 3 + YfE + i(a 2 + <S 2 )]s 2 + SEs 



(21) 



Equation ( |f9| ) describes the motion of a unit 'mass' with coordinate s = S z in a 'potential' U(s) so 
that its 'energy' E — |s 2 + U(s) is conserved. The true dimension of £, which is called e in Ref. |J, 
is (energy) 2 . In principle, we fix E by specifying s(0) and s(0). In practice [§, §2.2] there are only two 
independent constants of the three-dimensional motion of a spin, viz., the energy E and S 2 = er 2 , so that 
£ is bound to be a function of both of them: £ = ^(a 2 a 2 — E 2 ). This is most easily verified by using 
(pp|), computing E 2 , and realizing that E = —F(s) — aS x while s — —aS y . 

For E negative enough, U is a double-well potential, which is asymmetric in s for 5^0 and we 
have two disjoint closed orbits, as is brought out by Figs. 1 and 2. Figure 2 shows that the asymmetry 
develops much more slowly than the shift of the maximum; this latter is roughly at —8/2. The turning 
points a\ < b\ < 62 < 02 are solutions of the equation U(s) = £. Given E, classically allowed motion 
is between a 2 : and bi, with i = 1,2, so either on the left or on the right. The period of this motion and, 
thus, implicitly the attempt frequency is 



Ti{E) 



ds [2(£-U(s))} 



-1/2, 



(22) 



If £ is close enough to the bottom of the potential well, U(s) is nearly parabolic in the domain of integra- 
tion of d|) and Ti(E) w 2vr/ v /[/"(s;), where si and S2 are the locations of the minima. Because U(s) is 
independent of the spin quantum number S, the attempt frequency will also be (nearly) independent of 
S. 
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Once an orbit exists, any E is classically acceptable. Quantum-mcchanically, however, only 2S + 1 
energy eigenvalues survive. Determining the allowed eigenvalues to fair approximation and in closed form 
is what semiclassical quantization is (or should be) good for. 

Handling a single spin, we have only a single pair of canonically conjugate variables q and p, which are 
related through the Poisson bracket {q,p} = 1. Since in a Hamiltonian formalism q and p are handled on 
an equal footing there is no harm in interchanging them by putting q ncw := — p and p now := q so that the 
new variables have the same Poisson bracket {q n cw,Pncw} = 1 as the old ones. Instead of declaring q = S z 
and p = —(j) to be canonical coordinates [^|, 10 , with <f> as the azimuth, we now find it advantageous 
to put 

q = 4> and p = S z . (23) 
Semiclassical quantization is a condition on the action integral, 



pdq = <b S z d<j) = nh, n e 7L , (24) 

where the integral is to be taken over a classical, closed orbit, h is Planck's constant, and n is an integer; 
cf. Messiah |14{| . Along a closed orbit S z can often, certainly in the present case, be specified as a function 
of < 4> < 2ir; cf. Fig. 1. We write S z — <rcos8, with 9 as the polar angle, and arrive at the condition 
which we will use in the spin problem below, 

a(cos8) = ^- <j)d(f>cos9 — nh, n e 7L . (25) 

Of course cos 9 is to be given as a function of <p. The fact that n is a positive or negative integer, 
restricted to \n\ < S, is typical to spins. 

To verify that all this makes sense, we take the limit a — > so that the paraboloid — [jS z +5S z ] — aS x = 
E becomes very steep and, consequently, S z in ( §jj) is more or less constant as the spin tracks its orbit. 
Hence we find 2ttS z — nh, which is equivalent to saying S z — nh, as should be the case on the spectrum 
ofS 2 . 

As an application of the semiclassical quantization condition J25|) we study the Hamiltonian (j^) with 
vanishing S. The energy as given in polar coordinates, 

E = -7S2 - aS x = -7CT 2 cos 2 9 - acr sin cos <j> , (26) 

leads to a quadratic equation for y := x 2 with x — cos 9; this only happens when 6 = 0. Here we 
assume E < 0, which is typical to tunneling. In view of considerations to come in §4, we introduce the 
dimensionless parameters 



a ^, E , a cos 6 

— Q 2 = j and e=^^. 27 



Here a is supposed to be small; in the Mni2 case, a = 0.3. 

The quantity y obeys the equation y 2 - 2Q 2 (1 - 2e 2 Q 2 )y + (1 - 4e 2 )Q 4 = so that 

V± = Q 2 {(1 - 2e 2 Q 2 ) ± 2e [1 - Q 2 + e 2 Q 4 } 1 / 2 } > . (28) 

In view of y_(</>) = y+{4> + tt), h suffices to consider y + and, thus, 

x ± = ±Q{(1 - 2e 2 g 2 ) + 2e [1 - Q 2 + e 2 Q 4 ] 1/2 } 1/2 = ±Qf{Q, e) • (29) 

That is, we are left with a 'positive' branch, x+((j>), and a 'negative' one, x_ (0), symmetrically positioned 
with respect to the S x -S y plane. Both x+ and x- have to be inserted int o (|25|) . They correspond to 
n > and n < 0, respectively, and yield the same energy. Indeed, (^5| ) and fl29) imply 



g) 2 = Q 2 (/> 2 . (30) 
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As a matter of fact, / is a function of Q 2 so that the solution of ( |30| ) for E only depends on \n\. It then 
remains to calculate the quantum-mechanical splitting of the levels —n and n as indicated in §2. 

In the present case, semiclassical quantization is a straightforward integration giving up to second 
order in a (see Appendix A) 



(cos ( 



1 /" 27r 



±Q[1 



<Q- 



Q- 



rih 
<J 



n 
S 



(31) 



The dependence of x± (</>) on </> is a dependence upon cos </>. If desired, one can change variables through 
z := cxp(i0) and obtain a contour integral in the complex z plane. 

What we are after is the energy E as it appears in Q. First we solve ([3l]), a fourth order equation 
in Q. Recalling that it already contains an error term of the order of a 4 , it suffices to find Q = Qia) up 
to second order, which can be done by iteration. To this end we rewrite (^) in a form that is easy to 
iterate, 

Q= l -^ + ^(Q- 1 + Q- 3 )- (32) 

Iterating once, i.e., replacing Q in the right-hand side of (|32|) by H/S^ we obtain Q up to second order 



Q 



S 



a 

T 



(33) 



Because E = —jaQ, squaring ( p3| ) and dropping terms of order a we obtain the second-order expression 
for the frequencies (|n| < S), 



Q n = E n /h = -n 2 T - — 



1+ - 



(34) 



Here and in the next section we have dropped the superscript 'sc', and E n denotes a semiclassical 
eigenvalue. As announced, the degeneracy of E n and E_ n is not lifted. It has to be borne in mind that 
([34]) has been obtained under the assumption that we may drop everything beyond second order in a 2 , 
an assumption that may, but need not, hold. In Appendix A we derive an exact expression of (cos#) 
in terms of a power series in a. Using this power series one can obtain semiclassical energies up to any 
order. 

Equation ( p4|) being second-order in a, it is instructive to compare it with second order quantum- 
mechanical perturbation theory ||15[ , 



-n 2 T 



a 
8T 



n 2 + S(S+l) 



1/4 



(35) 



The correction to — n 2 T is slightly bigger in absolute value than in (|34|) but the agreement is excellent, 
except for, say \n\ < 2 and S < 4. For the ground state with |n| = 5, Eq. ( ]34| ) yields E±s = — ~/<t 2 ~ a 2 /47 
which agrees with the minimal classical energy, cf. ||, Eq. (2.6)]. Table I of Ref. ||] shows that, for 
a = 7 = ti = 1 and S > 8, the ground state has tts + S 2 T = —0.26, which is indeed near the predicted 
— 1/4. In fact, it is slightly less, as we would expect. Under the proviso S > 5 and |n| > 3, the deviation 
of the 'shift' E„ /ft + n 2 T from that given by numerically exact eigenvalues E n is less than 15%. 



4 Determining the Energies for Nonzero £ 

We proceed in analogy to the 6 = case. In the argument below we are looking for solutions of a fourth- 
order equation P4 (x) = in dependence upon a given combination a of the coefficients of the polynomial 
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P4. Instead of attempting to obtain an exact solution, which would not provide much insight, we rewrite 
Pi(x) = in the form of a fixed-point equation x = f(x,a), and, supposing the smallness of a, find the 
solution x — x(a) by iteration, up to a given order in a. We have already applied this procedure once, 
viz., to (||). 

The energy as given in polar coordinates and to be compared with (|26j) reads 

E = -jS 2 - SS Z - aS x = -'fa 2 cos 2 9 - 5a cos 9 - aa sin 9cos4> , (36) 
As before, we put x = cos 9 but do not get a quadratic equation in y := x 2 once 6^=0, 



ja 2 x 2 +Sax + aa\J 1 - x 2 cos + E = . (37) 
In agreement with (E7j) we now define the dimensionless quantities 



n 



<5 „ -E .o a cos< 



-r_ d=— Q 2 = ^ + rf 2 ^^7^ (38) 

270- 2717 70-^ 



which reduce to ( p7| ) whenever 5 = 0. Here too, i? < will be assumed. 
The two solutions x = x± of ( |37| ) obey the equation 

1/2 



d±Q(l-2eVl-a; 2 ) • (39) 



According to what has been outlined at the beginning of this section, we assume e < 1/2 so that the 
outer square root can be expanded. Then (|3^) can be solved by iteration to any order in a. The algebra 
has been relegated to Appendix B. After an integration with respect to one finds, to second order in a, 



If 2 * a 2 



(X EE 



2Wo ^ X ^= X ° + W [xh m 
where xq = —d ± Q and 

[x] 2 = T Q l±Qx + \{l - xl)] = T | (Q 2 + 1 - d 2 ) . (41) 

Unless stated otherwise, we will not repeat that henceforth we have to add a term 0(a 4 ) to all right-hand 
sides of the equations in this section. 

Semiclassical quantization means, in complete analogy to (|3l|), that (x) is to be equal to n[S with 
—S < n < S. Realizingthat n > corresponds to the upper and n < to the lower sign in (|4lJ) and 
±n = \n\, we combine @ with © and obtain 

2 

^ = {x)=±Q-dT^[l + Q- 2 (l d 2 )] (42) 



and thus 



^ = l 4 = Q Bgn(n) d-^[l + Q- 2 (l - d 2 )} . (43) 



We now rewrite this in a form that is apt to iteration, 



a 2 ^ n,. ,o„ „ ;/ 



Q = Qo + ^[l + Q- 2 (l-d 2 )} , Q = ^+sgn(n)d. (44) 



To obtain Q to second order in a, we iterate once 

,2 
h 



Q = Q» + ^ r li + Qo z (i-d 2 )} ■ (4--,) 
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For 8 = d = we recover (33). Squaring 
we obtain 

,2 



dropping terms of order a 4 and taking advantage of (B§ 



= E n /H — 



n r — rt<5 

sr 



iS/T 



0{a 4 TS 2 ) 



(46) 



{n + S/2T) 2 

Plainly, this is identical with (|34| ) for 5 = d = 0. On the other hand, second-order perturbation theory 



gives 



E^/h=-n 2 T- 



8T 



2 r 



5(5 + 1) 



nS/T 



(47) 



(n+ 5/2r) 2 - 1/4 

see also (JOm. The above expression agrees with Eq. ( |35|) when <5 vanishes, and with ( fi"6|) whenever n 
and S are sufficiently large and a < 1. The latter condition is quite reasonable since we have used 
quantum-mechanical perturbation theory and, thus, compared the 'perturbation' — aS x with the 'rest', 
viz., —7 S 2 — SS Z . 

Once we know the semiclassical energies E n , we can tune 8 so as to get semiclassical resonance 
E m = E n for some m 7^ n. Despite being obtained for a = 0, the resonance condition (^) remains valid 
for a as well. Both (46) and p^ ) depend on n through the expression n 2 + nS/T. Therefore (||) 
directly implies E m = E n for each couple {m, n} satisfying m + n = —k = —5/T. 

A closer inspection shows that (|[) implies degeneracy, to any order, of semiclassical eigenvalues. To 
see why, we return to Eq. (39). Let us start by averaging it, viz., 



n 
S 



-d±Q((l-2eVl~x 2 ) 1/2 } 



(48) 



We now add d on the right and on the left, multiply both sides by S, and square the result so as to find 

2 



n 2 +n8/r=-(S/2T) 2 + SQ((1 - 2e^l - x 2 ) 1 / 2 ) 



(49) 



The right-hand side does not show any n -dependence. It is obtained by taking x as a solution to the 
fourth-order equation (|39|), depending on E, a, 5, T, and (f>. After integration with respect to (f>, we are 
left with (f49|), an implicit equation for E. Solving it for E, the solution E = E n will depend on n through 
the combination n 2 + nS/T, as it shows up in ([49|). Hence E m — E n whenever m 2 + mS/T = n 2 + nS/T. 



5 Extensions 

The extension of the above results to half-integer spins is straightforward. All that we have to do is 
interpreting m and n as half- integers whenever they refer to eigenvalues of S z . In particular, m and n 
are half- integers in the resonance condition (§,§,||) and in the quantization condition ( p4] , ^5| ), and the 
semiclassical eigenvalues are also labelled by half-integers. At semiclassical resonance S/T is still an 
integer. 

The case where in the Hamiltonian S x is replaced by S l x with I > 1 brings nothing new, if I is an odd 
integer. For positive integer / semiclassical quantization is essentially unchanged. In the definiti on (|38|) 
of e, cos0 is replaced by cos' and in Eqs. (j39|), ( f4g| ) and (^) (%/l — x 2 ) 1 is substituted for Vl — x 2 . 
The upshot is that the semiclassical eigenvalues still depend on n through the expression n 2 + nS/T and 
therefore (|J) still implies E^ = En^ ■ If I is odd, spin tunneling splits the semiclassical degeneracy. 

For even I, however, the meaning of Eqs. (|^)-(|^) changes. Depending on the parity of fc, (|^) may 
refer to resonance or to true degeneracy. Let us recall that (^|) gives the condition for degeneracy when 
a = 0. In the case of integer spins, if k is odd then n — m — —k — 2m is odd for each m and n such 
that m + n = —k. As a consequence, the aS l x term does not split the degeneracy of all these level pairs. 
In the case of half-integer spins the same happens if k is even. Therefore E rn — E n holds for the true 
eigenvalues, even if not exactly at, but only very close to, the field S = kT at which E^ = En°\ 

An explanation of this generalized Kramers degeneracy can be found in §8 of |Q, but for the reader's 
convenience we repeat the argument here. We write the matrix of 7i (for I even) in the basis of the 
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iSz -eigenvectors |n). By noticing that (m\H\n) — whenever n — mis odd, we can permute the elements 
of the basis so as to transform the matrix into a block-diagonal form containing two blocks: the first is 
formed with n = — S, —S + 2, —S + 4, . . ., the second with n = —S + 1, —S + 3, —S + 5, . . .. Both can be 
diagonalized independently, so the von Neumann- Wigner argument [ |l2| about the avoided level crossing 
does not apply. 

For a degenerate energy level there is no tunneling. That is to say, a left-localized initial state 
remains approximately left-localized during time evolution. The reason is that in the corresponding two- 
dimensional eigensubspace of the Hamiltonian there exists an approximately left-localized state which is 
dominant in the eigenfunction expansion of the initial state. Therefore, if in an experimental situation 
the tranverse anisotropy fields are a combination of even powers only, measurements d la M, Ejj of the 
quantum decay of magnetization in an increasing longitudinal external magnetic field would ideally yield 
a fast decay (a step), only when S/T passes every other integer value. In other words, every second 
(expected) step would be missing. Unfortunately, an experiment in an applied magnetic field usually has 
a transverse component that induces steps at the remaining integer values. 

6 Discussion 

We have defined semiclassical resonance to be the degeneracy E m = E n for energy levels as they follow 
from semiclassical quantization. The ensuing level splitting is a purely quantum-mechanical phenomenon. 
Both can be handled straightforwardly by the WKB formalism || ^, |l(J . The resonance experiments as 
set up for spin tunneling in Mni 2 acetate and Feg magnetic crystals [|l6| are both elegant and quite 

promising since they provide a detailed test of the physics of spin quantum tunneling and the crossover 
from the regime of thermal activation to that of quantum behavior. The present paper's considerations 
on tunneling resonance are valid, if the temperature is low enough; say || Q below 1 K in the Mn 12 case. 

If so, one could apply semiclassical quantization (§3) to arbitrary numerical precision but the amount 
of insight thus obtained is fairly restricted. We have chosen a different way out and derived an explicit 
expression for the energy levels and, hence, for semiclassical resonance through an external field —5S Z 
under the assumption that a/(2TSQ 2 ) < 1 where Q 2 = (<5 2 /4 - jE)/(TS) 2 is dimensionless; cf. (||). 
For the ground state of Mni2 with S — 10, we have Q « 1 and the considerations apply. In fact, we 
have also derived the energy levels for S — under slightly less restrictive conditions and found that the 
shift of the ground state energy is well-predicted by — a 2 /4r. In some resonance experiments one would 
like to focus on states with m and n closer to zero since, in resonance, the tunneling frequency is much 
higher. As a consequence, the afore mentioned assumption may, but need not, hold. It mainly depends 
on Q 2 . If E w —n 2 TiT, then Q 2 i=s (n/S) 2 and it is the smaller the closer n is to 0. 

Whenever 6 is tuned to degeneracy of some energy levels of the Hamiltonian Hq = —jS 2 — SS Z , this 
degeneracy is not lifted by semiclassical quantization of the system with q^O. Only quantum mechanics 
does lift it, and only in high orders of a/TS. At a first sight rather surprisingly - but at a first sight only 
(§2) - the mathematics underlying tunneling for S ^ 0, when the Hamiltonian's symmetry of a rotation 
through 7r about the S x axis is broken, is identical to that of the case 5 — 0. That is to say, the very 
same formalism, whether WKB or otherwise, can be applied to both cases. 
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Appendix A: Semiclassical Quantization for 5 = 

We have to integrate x = cos 8 relative to <j>, i.e., average, the expression {. . .j 1 / 2 in j2S| ) and below, 

x ± = ±Q{1 - 2e 2 Q 2 + 26 [1 - Q 2 + e 2 Q 4 } 1/2 } 1/2 , (50) 
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where e = Q~ 2 a cos <f>. For 5 = 0, there are two disjoint classical orbits once E < —aa. Exactly in this 
energy range one finds {...}> so that taking the square root in ([30|) is well-defined. For the moment 
we suppose that a and thus e is small enough to perform a series expansion of the outer square root. It 
is of the form 

(2n-3)!! 



{I- 2V) 1 ' 2 = 1-J2b n y n where b n 



(51) 



The series converges for \y\ < 1/2 while b\ = 1, b% = 63 = 1/2. Introducing the short-hand r 
[1 — Q 2 + (eQ 2 ) 2 ] 1 / 2 we see that we have to average 



{1 - 2e(eQ 2 - r)} 1 ' 2 = 1 - £ b n e n (eQ 2 - r) n . 



n>l 



(52) 



We now concentrate on the series, insert the definition of e and apply the binomial theorem to (|52|) so as 
to find 

n 



b n Q- 2n a m+n (-iy 



n>l rn—0 



(53) 



At this point we observe that r depends on <f) through cos 2 </>. Hence, (g(r) cos e 4>) vanishes for I odd, 
whatever the function g may be. Thus, when averaging (|5^), only terms with I ~ m + n even will survive. 
For these terms n — m = (n + m) — 2m is also even and we can apply the binomial theorem to r n ~ m as 
well. Furthermore, averaging a cosine is equivalent to a contour integration over the unit circle Ci in the 
complex z plane, 

(54) 



2 i" dz , -ut 
(k —(z + z 



Finally, we are left with 



(cos0) = ±Q 



2ni 



n (n—m)/2 



£/2 



!>lm=0 (=0 



(n — m)\ ( 21 + m + n 



I J\l + Um + ri) 



±Q 



21: 



k=l 



(55) 



where it is understood that m + n is even. If it were not because of the third binomial coefficient in 
(|S5|) we could directly resum the series. In the style of Borel |i~7|, Q, however, we could, by noting 

|~0[ for the Gamma 

for its inverse, 



n\ = T(n + 1) and taking an Eulerian integral representation of the second kind |19 
function in conjunction with Hankel's contour integral 



1 



T(t) 



1 

2^ 



(-*)- 



(56) 



where C2 is a counter-clockwise infinite contour around the positive real axis. We will not pursue this 
idea here. Collecting terms of order a 2 we find C\ = Q~ 2 + Q~ 4 and arrive at (|3l|), as advertised. 

How good is all this? Since we interchange the series ( |5l|) and averaging, viz., § dz . . ., we have to 
prove uniform convergence of the series and thus, in view of (|5l| ), uniform boundedness of y = e(eQ 2 — r) 
as cos 4> varies between 1 and — 1 . The upshot appears in (pa). Plainly, the series in (|55|) is convergent 
for a small enough so that, for a given a, Q and, thus, the energy has to stay away from zero. The series 
in (55) certainly diverges at energies for which the two classical orbits merge into a single one. One can 
handle divergent series through, e.g., Borel summation |l7], |l8| as above, but this is not the topic of the 
present paper. 
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Appendix B: Quantizing to Second Order in e 

In this appendix we consider e as a 'small' parameter and sketch how one can perform semiclassical 
quantization up to second order in e. As we will see, the method can be generalized and the considerations 
of Appendix A concerning the validity of the approximations involved apply here as well. 



We denote the two branches x± of (39), viz 



x = -d±Q[l- 2e\/l~x 2 } 1 ' 2 , (57) 

simply by x. Keeping in mind that x = cos 8 depends on <f>, we fix <f> for the moment and study x's 
dependence upon e = aQ~ 2 cos0. In view of ( |57j ) it is clear that x — x(e) is analytic in a neighborhood 
of e = 0. Defining the nth partial 'Taylor sum' of x(e) to be 

^E^r^EM^, (58) 

k=0 ' k=0 

with Xq = x(0) — — d ± Q, using (|5l|), and expanding the outer square root of (p7|), we can rewrite (|57|) 
in the form 

n>l 

Up to and including terms of order two we then find 

x = x T Q[e(l - x 2 ) 1 ' 2 + e 2 (l - x 2 )/2] + 0(e 3 ) , (60) 

whence 

x 1 ^xoTQ(l-x 2 ) 1/2 e = x + [x} 1 e (61) 

because y/l — x 2 = ^/l — x 2 , + 0(e). 

In order to obtain x 2 , we observe that (pOI) implies 



[x} 2 = TQ{[Vl - x 2 V + [1 - x\ /2} . (62) 



In analogy to (|5g)) we have used the notation [/]„ = f^(0)/n\. 
As for the entries of ([62]), a simple calculation gives 



[Vl-a; 2 ]! 



dvT 



de 



= -x (l-^)-VV(0) 
-z (l-Zo)" 1/2 Mi =±Q^o ■ (63) 



Since [1 — x 2 ]q = 1 — x\ we end up with 



:r 2 = xi + [x] 2 e 2 = xqT Qe\J 1 - x 2 - Qe 2 [Qx ± (1 - x 2 )/2] (64) 

It is £2 that is to be averaged with respect to </>. Inserting (e) = and (e 2 ) = a 2 /2Q A into (x 2 ) we obtain 
Eq. ©. 
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